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employing a microscopic point of view. A comparison is also made between our 
results and those of other many-body techniques. 
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1 Introduction 

Pulsars are rapidly rotating neutron stars with strong surface magnetic fields in the range 
of 10^^ — 10^^ Gauss [1, 2, 3]. The physical origin of this magnetic field remains an open 
problem and there is still no general consensus regarding the mechanism to generate such 
strong magnetic fields in a neutron star. There exist several possibilities of the generation 
of the magnetic field in a neutron star, from the nuclear physics point of view, however, 
one of the most interesting and stimulating mechanisms which have been suggested is the 
possible existence of a phase transition to a ferromagnetic state at densities corresponding 
to the theoretically stable neutron stars and, therefore, of a ferromagnetic core in the 
liquid interior of such compact objects. Such a possibility has been studied by several 
authors using different theoretical approaches [4-25] , but the results are still contradictory. 
Whereas some calculations, like for instance the ones based on Skyrmelike interactions 
predict the transition to occur at densities in the range (1 — 4)po (po = 0.16/m~^), others, 
like recent Monte Carlo [20] and Brueckner-Hartree-Fock calculations [21-23] using modern 
two- and three-body realistic interactions exclude such a transition, at least up to densities 
around five times po- This transition could have important consequences for the evolution 
of a protoneutron star, in particular for the spin correlations in the medium which do 
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strongly affect the neutrino cross section and the neutrino mean free path inside the star 
[26]. 

In recent years, we have computed the equation of state of symmetrical and asymmet- 
rical nuclear matter and some of their properties such as symmetry energy, pressure, etc. 
[27-30] and properties of spin polarized liquid ^He [31] using the lowest order constrained 
variational (LOCV) approach. The LOCV method which was developed several years ago 
is a useful tool for the determination of the properties of neutron, nuclear and asymmet- 
ric nuclear matter at zero and finite temperature [27-39]. The LOCV method is a fully 
self-consistent formalism and it does not bring any free parameters into calculation. It 
employs a normalization constraint to keep the higher order term as small as possible 
[34, 27]. The functional minimization procedure represents an enormous computational 
simplification over unconstrained methods that attempt to go beyond lowest order. 

In the present work, we compute the polarized neutron matter properties using the 
LOCV method with the AVis potential [40] employing microscopic calculations where we 
treat explicitly the spin projection in the many-body wave functions. 

2 Basic Theory 
2.1 LOCV Formalism 

We consider a trial many-body wave function of the form 

= Fcf>, (1) 



3 



where is the uncorrelated ground state wave function (simply the Slater determinant 
of plane waves) of independent neutron and F — F{1 ■ ■ ■ N) is an appropriate N-body 
correlation operator which can be replaced by a Jastrow form i.e., 

F = Sl[f{zj), (2) 

i>j 

in which S is a symmetrizing operator. We consider a cluster expansion of the energy 
functional up to the two-body term, 

mf\) = i^ = E, + Er (3) 



The one-body term Ei for a polarized neutron matter can be written as Fermi momentum 
functional = (GttV^'^)^): 



Labels 1 and 2 are used instead of spin up and spin down neutrons, respectively, and 
p — p^^^ -I- p(^) is the total neutron matter density. The two-body energy Ei is 

^2 = ^E(^j>(12)|ij-ji), (5) 

where 

1/(12) = -£[/(12), [V?2, /(12)]] + /(12)y(12)/(12), /(12) and ^(12) are the two-body 
correlation and potential. For the two-body correlation function, /(12), we consider the 
following form [27, 28]: 
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fe=l 



/(12) = ^/«(12)0(^)(12), (6) 



where, the operators 0^'^\l2) are given by 



(7) 



and 5*12 is the tensor operator. 

After doing some algebra we find the following equation for the two-body energy: 



_ 2 (h'\ (2J + 1) 

2m 



dr 



.'2 



2772 

+7^({^c +Va + Vr + + (-6A; + U){Vu + V^t) - (A; - + V^) 

+[Vt + Kr. + (-6A; + U)VtT\[2 + 2Vr,]}a^^^\k^r) 

+[Vn + V^2. + Vi2r + 1^2.r]cL'=)'(A;/r) + [\^,2 + V;,2r]4'^'(A;/r))/(^)'] 



(8) 



where a = {J, L, S, Sz} and the coefficient a^^^^, etc., are defined as 



(9) 



(10) 



(11) 



(12) 
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c(f)^(x) = [773/7-1,5. (2:) + ^^3/7+1,5.(2:)], 



d^^\x) = x'[C2lj-i,sA^) + A2/.+l,5.(^)], 



d^fix) = Ai3lj-1,SAX) + A3/J+1,5.(X)], 



^ 2J+1' ^ 2J + 1' 2J+1 ' 



+ "^ = ^jT1-' = 2JTI ' 



_ J(J^ + 2J+1) _ J(J^ + J + 2) 
~ 2JTi ' ~ 2JTi ^ 



+ 2 + 2 J + 1 _ J(J^ + J + 4) 
~ 2J+1 ' '^^ ~ 2J+1 ^ 



_ J(J^ + J+1) _ + 2 + 5 J + 4 



/j,5.(x)= JdqPsMUjM- 
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In the above equation Jj{x) is the Bessel's function and, Ps^iQ) is defined as follows, 

for I - /c^^^l <q<l + k''/'\, 

PsM) = l'^mMkF'\k''F'), (25) 

for q < l\k'^'' -k^/^\ and 

PsAq)-0, (26) 

for g > I j/cp ^ + k'p^l, where cr^i or az2 — |) ~| for spin up and spin down, respectively. 
Now, we can minimize the two-body energy Eq.(8), with respect to the variations in 
the function fa^^^ but subject to the normalization constraint [28], 



(27) 



where in the case of spin polarized neutron matter the function hs^{r) is defined as 

-1/2 

hsXr) = 1-9 



k\ 



; -5, = ±1 

; 5, = 



(28) 



Prom the minimization of the two-body cluster energy, we get a set of coupled and un- 
coupled differential equations which are the same as presented in Ref.[28]. 
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2.2 Magnetic Susceptibility 

The magnetic susceptibility, x, which characterizes the response of a system to the mag- 
netic field, is defined by 

^=M„.o- '''' 
where M is the magnetization of system per unit volume and H is the magnetic field. By 
some simplification, the magnetic susceptibility can be written as 

where /i is the magnetic moment of the neutron and 5 is the spin polarization parameter 
which is defined as 

_ p(2) 

S^^ (31) 

P 

Usually, one is interested in calculating the ratio of x to the magnetic susceptibility for 
a degenerate free Fermi gas (xf). Xf can be straightforwardly obtained from Eq. (30), 
using the total energy per particle of free Fermi gas, 

where kp — (37r^p)^/^ is Fermi momentum. After a little algebra one finds 

X 2 Ef 



where Ep — h^kp/2m is the Fermi energy. 



(33) 
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3 Results 

In Fig. 1, we have shown the energy per particle for various values of spin polarization of 
the neutron matter as a function of density. As can be seen from this figure, the energy of 
neutron matter becomes repulsive by increasing the polarization for all relevant densities. 
According to this result, the spontaneous phase transition to a ferromagnetic state in the 
neutron matter does not occur. If such a transition existed a crossing of the energies of 
different polarizations would be observed at some density, indicating that the ground state 
of the system would be ferromagnetic from that density on. As is shown in Fig. 1, there 
is no crossing point. On the contrary, it becomes less favorable as the density increases. 
For the energy of neutron matter, we have also made a comparison between our results 
and the results of other many-body methods with the AVis potential as shown in Fig. 2. 
The BGLS calculations are based on the Brueckner-Hartree-Fock approximation both for 
continuous choice (BHFC) and standard choice (BHFG) [41]. The APR results have been 
obtained using the variational chain summation (VCS) method [42] and the EHMMP 
calculations have been carried out using the lowest order Brueckner (LOB) technique 
[43]. We see that our results are in agreement with those of others, specially with the 
APR and EHMMP calculations. 

For the neutron matter, we have also considered the dependence of energy to the spin 
polarization 5. Let us examine this dependency in quadratic spin polarization form for 
different densities as shown in Fig. 3. As can be seen from this figure, the energy per 
particle increases as the polarization increases and the minimum value of energy occurs 
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ai 5 — for all densities. This indicates that the ground state of neutron matter is 
paramagnetic. In Fig. 3, the results of ZLS calculations using the Brueckner-Hartree- 
Fock theory with the AVis potential [25] are also given for comparison. There is an 
agreement between our results and those of ZLS, specially at low densities. From Fig. 3, 
it is also seen that the variation of the energy of neutron matter versus 5^ is nearly linear. 
Therefore, one can characterize this dependency in the following analytical form 

E{p,S) = E{p,0) + a{p)S'- (34) 

The density dependent parameter a{p) can be interpreted as the measure of the energy 
required to produce a net spin alignment in the direction of the magnetic field, and its 
value can be determined as the slope of each line in Fig. 3, for the corresponding density, 

dE{p,S) 

<P) = • (35) 

In Fig. 4, the parameter a{p) is shown as a function of the density and as can be seen 
the value of this parameter increase by increasing the density. In turn this indicate the 
energy which require to align spin at same direction increases. An conclusion can be 
inferred again from this result is that a phase transition to a ferromagnetic state is not to 
be expected from our calculation. The parameter a{p) obtained by ZLS [25] is also shown 
in Fig. 4, for comparison. 

In Fig. 5, we have plotted the ratio x/Xf versus density. As can be seen from Fig. 
5, this ratio changes continuously for all densities. Therefore, the ferromagnetic phase 
transition does not occur. For comparison, we have also shown the results of ZLS [25] in 
this figure. 
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The equation of state of polarized neutron matter, P{p,S), can be simply obtained 



using 



(36) 

In Fig. 6, we have presented the pressure of neutron matter as a function of density p 
at different polarizations. We see that the equation of state becomes stiffer by increasing 
the polarization. We also see that with increasing density, the difference between the 
equations of state at different polarization becomes more appreciable. In order to check 
the causality condition for our equations of state, we have calculated the velocity of sound, 
Vs, as shown in Fig. 7. It is seen that the velocity of sound increases with both increasing 
polarization and density, but it is always less than the velocity of light in the vacuum (c) . 
Therefore, all calculated equations of state obey the causality condition. 

As it is known, the Landau parameter, Gq, describes the spin density fluctuation in the 
effective interaction. Gq is simply related to the magnetic susceptibility by the relation 

^ = (37) 

where m* is the effective mass. A magnetic instability would require Go < —1- Our 
results for the Landau parameter have been presented in the Fig. 8. It is seen that the 
value of Go is always positive and monotonically increasing up to highest density and does 
not show any magnetic instability for the neutron matter. In Fig. 8, the results of ZLS 
calculations [25] are also given for comparison. 
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4 Summary and Conclusions 

The properties of neutron matter is of primary importance in the study of neutron star, 
and in particular, strongly magnetized ones (i.e. pulsars). It is therefore important to 
calculate the properties of polarized neutron matter, using an efficient and sufficiently ac- 
curate method. We have recently computed various properties of the neutron matter using 
the lowest order constrained variational (LOCV) scheme. In order to make our results 
more general, we used this method for the polarized neutron matter. Energy per particle 
for various values of spin polarization of the neutron matter was computed as function 
of density, and shown to become repulsive as a result of increasing the polarization. In 
addition, we considered the dependence of energy of neutron matter to the spin polariza- 
tion, and found it to increase with the spin polarization for all densities. This dependence 
was represented by a quadratic formula with the coefficient of the quadratic term a{p) 
determined as a function of the density. This parameter, too, was shown to increase 
monotonically with density. Magnetic susceptibility, which characterizes the response of 
the system to the magnetic field was calculated for the system under consideration. We 
have also computed the equation of state of neutron matter at different polarizations. 
Our results for higher values of polarization show a stiff equation of state. The velocity 
of sound was computed to check the causality condition of equation of state and it was 
shown that it is always lower than the velocity of light in vacuum. We have also inves- 
tigated the Landau parameter Gq which shows that the value of Gq is always positive 
and monotonically increasing up to high densities. Finally, our results showed no phase 
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transition to ferromagnetic state. We have also compared the results of our calculations 
for the properties of neutron matter with the other calculations. 
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Figure 1: The energy per particle versus density (p) for different values of the spin polar- 
ization (6) of the neutron matter. 
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— ■ — Ours (LOCV) 
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Figure 2: Comparison between our results for the energy per particle of neutron matter 
and those of BGLS [41], APR [42] and EHMMP [43] calculations with the AVis potential. 
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Figure 3: Our results (full curves) for the energy difference of polarized and unpolarized 
cases versus quadratic spin polarization (5) for different values of the density(p) of the 
neutron matter. The results of ZLS [25] (dashed curves) are also presented for comparison. 
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Figure 4: Our results (full curve) for the parameter a{p) as a function of the density(p). 
The results of ZLS [25] (dashed curve) are also given for comparison. 




Figure 5: As Fig. 4, but for the magnetic susceptibility (x/Xf)- 
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Figure 6: The equation of state of neutron matter for different values of the spin polar- 
ization (S). 
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Figure 7: The velocity of sound in the unit of c versus density (p) for different values of 
the spin polarization (5). 
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Figure 8: As Fig. 4, but for the Landau parameter, Gq 
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